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Abstract 

We discuss tensor metric perturbations in a magnetic field around 
the homogeneous Juttner equilibrium of massless particles in an expand¬ 
ing universe. We solve the Liouville equation and derive the energy- 
momentum tensor up to linear terms in the metric and in the magnetic 
field. The term linear in the magnetic field is different from zero if the 
total charge of the primordial plasma is non-zero. We obtain an analytic 
formula for temperature fluctuations treating the tensor metric perturba¬ 
tions and the magnetic field as independent random variables. Assuming 
a cutoff on large momenta of the magnetic spectral function we show that 
the presence of the magnetic field can discriminate only low multipoles 
in the multipole expansion of temperature fluctuations. In such a case 
the term linear in the magnetic field may be more important than the 
quadratic one (corresponding to the fluctuations of the pure magnetic 
field). 


1 Introduction 

The magnetic field is ubiquitous in the universe. In particular, the CMB results 
from quantum thermal fluctuations of the electromagnetic field. It is present in 
the standard model. Fluctuations of the magnetic field may be expected in any 
model of the early universe. The non-trivial question concerns the appearance 
of a macroscopic magnetic field. There are various mechanisms which can be 
responsible for this phenomenon [1] [2]. There is no convincing argument for any 
of them. Let us mention the one which assumes a non-zero total charge of the 
primordial plasma[3]; the assumption relevant for this paper. 

The fluctuations lead to a diffusion of particle motions [4] and to a random 
rhs of the Einstein equations resulting from the energy-momentum. The energy- 
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momentum contains a contribution of the free electromagnetic field. This term 
has been studied in [5] [6] [7] [8]. When the magnetic field is Gaussian then the 
noise coming from the energy-momentum of the free electromagnetic field being 
quadratic is non-Gaussian. Non-Gaussian effects have not been discovered in 
CMB yet. This may be so because the magnetic field is weak and the quadratic 
terms are small. We point out in this paper that a particle interaction with the 
magnetic field leads to a contribution to the energy-momentum of the primordial 
plasma which is linear in the magnetic field. This happens if the total charge of 
the primordial plasma is non-zero. In such a case the impact of the linear term 
may be stronger than the one coming from the energy-momentum quadratic in 
the magnetic field. The strength of this term depends on the charge. There 
are strong bounds limiting the charge of the universe [9] [10]. We calculate a 
variation of the metric corresponding to the term depending on the magnetic 
held. We obtain an analytic formula for the temperature fluctuations result¬ 
ing from primordial fluctuations of the metric and of the magnetic held. The 
temperature fluctuations contain an information on structure formation. The 
impact of the primordial magnetic held on structure formation (and tempera¬ 
ture fluctuations) is usually ignored. There are however some arguments ( see, 
e.g. [11],p.575) indicating that the magnetic held should be taken into account 
in the studies of structure formation. 

The plan of this paper is the following. In sec.2 we find a perturbative so¬ 
lution of the Liouville-Vlasov equation describing a stream of particles in an 
inhomogeneous expanding metric and in the magnetic held. We are interested 
in the ultrarelativistic limit when all the particles are massless. As a zero order 
solution we choose the Jiittner distribution [12] with a time dependent temper¬ 
ature. In sec.3 we discuss Einstein equations with the energy-momentum on 
the rhs which is determined by the solution of the Liouville-Vlasov equation. A 
perturbative solution of Einstein equations determines a variation of the metric 
in the magnetic held. Fluctuations of the temperature are calculated in sec.4. 
Temperature fluctuations are expanded in Legendre polynomials (multipole ex¬ 
pansion). We study a dependence of the expansion coefficients on the spectral 
function of the stochastic primordial magnetic held. In the Appendix we discuss 
some technical aspects of the estimates on the spectral function of fluctuations 
of the magnetic held. 

2 Liouville-Vlasov equation 

In this section we solve perturbatively the Liouville-Vlasov equation describing 
a distribution of classical trajectories (see [13] [14] for its application in general 
relativity). We decompose 

9fiu — V h[4i/ : (1) 
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where h describes homogenous metric in the conformal time and 

ds 2 = g^dx^dx 1 ' = a 2 (dt 2 — dx 2 — 7 , jdx l dx^). ( 2 ) 

In eq.(2) we assume that the tensor perturbations 7 ^ are transverse and trace¬ 
less. We write the Liouville equation in the form 

(p**^ - T^py'dkWe = 8T k v p»p v d k Sl e + eF^g^djQe, (3) 

where e is the electric charge, T£ are Christoffel symbols, d* = and x = 
(t,x) (boldface letters denoting the three vectors), d k = denotes derivatives 
over momenta. 

For massless particles ( m = 0) and in the homogeneous metric (h^ = 
0) any function /(a 2 |p|) is a solution of eq.(3)(5T = 0)[15]. Because of the 
thermodynamic interpretation we choose the equilibrium distribution f Ie [ 12 ] 
as a starting point of the perturbation 

= 5(2tt ) _3 ( exp(a 2 /?|p| - pp(e)) + q) (4) 


with 

P 2 = 'Y^P’V 3 - 
3 

g is the number of particle’s degrees of freedom (we set g = 1 from now on), 
ft = ^ is the inverse temperature and p(e) is the chemical potential for particles 
of the type e. In eq.(4) q = 1 for fermions, q = — 1 for bosons and q = 0 for the 
classical Jiittner distribution. The physical momentum is a|p| and the physical 
temperature (pa)- 1 . 

We write 

fl e = n^(l + Xe) (5) 

and look for a perturbative solution \ of eq.(3) in an inhomogeneous metric (2). 
Then, 


dtX + n k d k x ~ 2'Hp k dkX + e |p| 1 F ju h UIJ ,p^d j x fr) 

= -a 2 Pf\p\(dajkn°n k + djji k n l n k ) + ea 2 Pf\p\~ 2 F : > v h VIJ ,p IJ 'p j , n>l 


where n k = p k \p\ 1 ,H = a 1 ^| and / = fV(a 2 /3|p|) where 

Cl'(x) = exp(:r — p(e)j3) ^ exp(x — g(e)(3) + q ^ 

/ = 1 for the Jiittner distribution (relativistic equilibrium distribution neglect¬ 
ing the quantum statistics). 
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We set / — 1 and look for a solution of eq.( 6 ) which is of the first order in 
momentum 

Xe = a 2 /3v\p\ +/3r e . (7) 

We assume an infinite conductivity of the primordial plasma. Then, the electric 
field is zero [16]. Inserting (7) in eq.( 6 ) we obtain equations for v and r 

d t v + n k d%v = -da jk nPn k - \n^d^ lk n l n k , (8) 

d t r e + n k d%r e = e^ km e jkl n j n m B l = a e , (9) 

where we wrote 

F jk = F kl Bi(t) 

with 

B(t, x) = a~ 2 B(x) 

This time-dependence of the magnetic field follows from Maxwell equations in 
an expanding universe [17]. We introduce 

B = jk n J n k ( 10 ) 

Then 

d t Q + n k d%Q = -\da 3 krin k = U (11) 

0 as discussed in our earlier paper [18] has the meaning of the temperature 
variation. The solution of eq.(ll) reads 

O t (x) = f* ds7l(s,x-(t-s)n), (12) 

The Fourier transform of the solution of eq.(9) is 

f e (k, t) = fg dsa e (k, n, s) exp(— ikn(t — s)) (13) 

We calculate (in the conformal time) the energy-momentum tensor till the 
first order in the metric and the magnetic field perturbation ( we preserve the 
quadratic term Tem of the free electromagnetic field for a later comparison). 
For this purpose we sum the densities of +1 and —1 particles [19] 

T jl = (2TT)- 3 ^gfd P Po yy(fl+ + fi_) +T% m 

= J dnn J n l (^l + 4v)(fl ( ^ ) + + t em ( 14 ) 

= ^ + ^ + 4 

where 

t em = - V l B k B k ) 

Here, Tq is the energy-momentum of the solution (4) of the Liouvillc equation 
on the homogeneous space-time, 5T denotes the terms linear in the metric 7 
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and in the magnetic field B. In eq.(14) v should still be expressed by 0 from 
eq.( 10 ) or determined from eq.( 8 ). 

We can calculate the tensor (14) taking multiple derivatives 

f d d 

/ dmexp(— ikm(f — s))m n ....m 3r = i r 7 —-. —si n (<?), (15) 

J oq n oqir 

where after the calculation of derivatives we should set 

q = (t — s)k. 

3 Einstein equations in a magnetic field 

We write the spatial part of Einstein equations for traceless transverse metric 
perturbations in the form (G® is the Einstein tensor, G is the Newton constant) 


6 G 3 = - 2 a 2 (5^7 ij - 3A7 ij + / Hdt'Yij^j 
= 8 na 2 GP. irM ( ST kl + ( T EM ) lk ), 


(16) 


where 

Pij\kl — AjfcA/j T Aj/Afcj Aij'A/fc, 

A jk — 5j k Oj d k A 

and A is the three-dimensional Laplacian. The matter energy-momentum 5T 
on the rhs is linearly dependent on the magnetic field whereas the free electro¬ 
magnetic energy-momentum Tem is quadratic in the magnetic field. 

We write 

7jfc = a _1 7jfc- (17) 

Then, eq.(16) takes the form 

Qljk = dtjjk - (A + a -1 G> 2 a)7jfc = 8 tt a 5 G((6T TT ) jk + (Tf^) ife ), (18) 

where the rhs still depends on the metric. For a general a(t) it is not simple to 
solve eq.(18). Let Q~ 1 {k\t,s) be the kernel of the inverse of Q. Then, we can 
solve eq.(18) by iteration (till the first order in G) 

jjk = rW + 8 ^ f dsQ-\k-,t,s)a 5 STj k T {s,a~ x 7 ^), (19) 

Jo 

where 7 9rav is the solution of the homogeneous equation ( at G = 0). 

In the radiation era d^a = 0. Then, 

C7 —1 (At; t, s ) = k~ x sin(fc(i — s)). (20) 

In another limit, if ka « 8%a then the dependence of Q~ x on k can be neglected. 


5 




Returning to eq.(18) we perforin some integrals over n (using eq.(15)) and 
write it in the form 

dtlij - Ajij + 2'Hdt'Yij = 87 t Ga 4 P ij;k iT^ l M , , 

+8nG(2'7r)~ 3 96na~ 2 p~ 4 ‘P i j-ki(^ki +t lk ), 

where 

t jl = j dnn j n l (& + 3(N+r+ + !V_r_)) = 9 jl + r jl . (22) 

0 is determined from eqs.(ll)-(12) and r is defined in (13) (6 and r do not 
depend on a). For the equilibrium distribution fl p we have assumed the ap¬ 
proximate formula (justified for high energies) 

= ^exp(-/la 2 |p|), (23) 


where from eq.(4) at q —> 0 


Nf, = exp(/z(e)/3) 

In order to write down Einstein equations explicitly we calculate the part of 
the plasma energy-momentum ( 22 ) which is linear in the magnetic field 

r jfc (k, t) = (N + — 7V_) f* dsdmm- J m fc cr(k, m, s) exp(— *km(t — s)) 

= (N + — N-) f* dsdmdqm : ’m k 'yii(q, s)e rlp B p ( k — q) exp(— ikm(t — s))m l m r . 

(24) 

m denotes the directional vector of propagation (which we denoted by n in 
eq.(9)). So, the part of the metric perturbation coming from the magnetic field 
is 


dtlij ~ Ayij + ‘2'Hdt'Yij = 8TrG(2n) 3 96na 2 f3 4 Pij-kir lk . (25) 
For the remaining part of eqs.(21)-(22) we have 

6T jk = (27r) _3 967ra _6 /3 _4 t jfc 

This part of metric fluctuations is discussed in many text-books [20] [21] (we 
have calculated it for a diffusive matter in [18]). 

4 Temperature fluctuations 

There will be temperature fluctuations caused by the density fluctuations (scalar 
perturbations), gravitational waves (quantum metric fluctuations) as well as 
fluctuations of the primordial magnetic fields. The solution for the temperature 
fluctuations is expressed by 7 (eq.( 12 )) 

B(£,n) = -nVi f* dsd s jji(s,x - (t- s)n), 
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where 7 is determined by P l . Now, 


(@(t, n)0(i, n')} = 

= \{2i t)~ 3 f* ds fg ds' f dkdk'n l nPn' r n' p (26) 

d s d s ' {'j j i{s,k)'y rp (s',k')) exp(-i(t - s)nk + i(< - s')n'k), 

In eqs.(25)-(26) we wish to calculate the part of fluctuations coming from the 
magnetic field. It is determined from the solution of Einstein equations (16). 
From the perturbative solution (19) the fluctuations coming from the magnetic 
field are 

( 7 jJ (s,k) 7 ob (s',k')) = ^ 87 rG( 27 r) - 3 967 ra -2 /? -4 ^ (a(s)a(s ')) -1 fg fg drdr' 
a(r)a(T')t/ _ 1 (fc, t)!/” 1 ) k ', s', T')Pji. p qP ab . rnn (r pq (a~ 1 7 , k, r)r m „(a _1 7 . k', t')), 

where r is expressed by the metric and by the magnetic field in eq.(24). 

We have two random fields in the solution 0: 7 and B. We assume that B 
is a Gaussian random field with the covariance 

(B^B^k')) = A iJ -(k)5(k + k , )P B (k), (27) 

where 

Aji(k) = Sji - kjkikT 2 . 

B( k) is time-independent as explained below eq.(9) ( we denote a function and 
its Fourier transform by the same letter; the meaning should follow from the 
context), 7 is an independent random field with the covariance 

( 7 j;(s,k) 7 afe (s , ,k')) = Pji. ab 5(k + k')P^(k; s, s'). 

We wish to calculate the correction to the temperature fluctuations coming 
from the interaction with the magnetic field. Let us define 

Peri q, S, s')S( q + q') = (cr(s, q)cr(s', q')}, (28) 

where 

er(k, m, s) = f dp 7 ifc(p, s)e rla B a (k - p )m k m r (29) 

Then, 

(cr(k, m, s)cr(k', m', s')) = / dpdp'( 7 ii(p, s)e rla B a (k - p)m l m r 
7;'j'(p , i s')e rla ' B a >{k' — p')m' r ' m' i ) 

= <5(k + k')A aa f (k)P H .jf v (k)e rla m l m r e a ' r ' 1 'rn ,r 'm H ' f dpP 7 (p; s, s')P B (k - p) 

= d(k + k')tu(k, m, m') / dpP 7 (p; s, s')Pb{ k - p), 

(30) 

where 

w(k, m, m') = A aa f{k)Pi i .if i f{k)e rla m l m r e a r 1 m' r m n . 
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The spectral function defining the fluctuations of eq.(26) is determined by 


P CT (k;s, s') 


dpPj(p', s, s')P B (k — p). 


( 31 ) 


We consider a simplified version of the graviton correlation function 


P 7 (fc; s, s') = f(s, s')k~ a 


(32) 


with 0 < a < 3 and 


P B {k) = b 2 k° exp(-y) 


(33) 


(usually with — 3 < a < 2) as a model for the spectral function of a primordial 
magnetic field with a Debye frequency cutoff for k 2 > A [17] [22]. We discuss the 
spectral function P a in the Appendix. We show that for a small k the spectral 
function P a (k;s,s') tends to Kf(s,s') (with a certain constant K) whereas 
for a large k it behaves like the graviton spectral function,i.e.,as k~ a . Hence, 
the magnetic field substantially changes the powerlike behaviour of the power 
spectrum for a small k but it does not change the leading behaviour for a large k. 
For comparison the contribution of the pure electromagnetic energy-momentum 
tensor to the spectral function of the temperature fluctuations is determined by 


P 2 u(k) = JdpP B (p)P B (k-p) (34) 

We show in the Appendix that for a > — 1 it tends to a constant for a small k 
and decays exponentially for k 2 > A. So, it does not contribute to high mul¬ 
tipoles. For low multipoles its contribution behaves as b 4 (the fourth power of 
the strength of the magnetic field) whereas the contribution to the tempera¬ 
ture fluctuations (24) of the term linear in the magnetic field is proportional to 
b 2 (N + — N_) 2 . Now, using (24),(26) and (30)-(31) 


(d t jji(k, = (8 ttG(2tt) 3 96na 2 /3 <5(k + \d)d t d t > fo ds f* ds ' 

(a(t)a(t'))~ 1 Q~ 1 (k-,t, s)Q~ 1 (k'-,t',s')a(s)a(s') f* dr / Q s dP 

Wji. jU '{s - r, r, s' - P, r'; k)P a (k; r, t 1 ) = 6(k + k')Pji. fll {k)F(t,t';k), 

(35) 

where 

^ r, r, s' - P, P; k) = / dmdm'Pji. a b(k)P j 'i>. a >^(k)m a m b m ,a 'rn ,b ' 
exp(ikrn(,s — r) — ikm'(s' — r'))w(k, m, m') = 

f dmdm'Pji. a b(k)Pj>i'. a 'b'(k)m a m b m' a m ,b exp(ikm(s — r) — ikm'(s' — P)) 
A 0 o' {k)Pu-i’i' (k)e rla m z m r e arl ' m’ r ' m n '. 

(36) 


The m integrals can be evaluated from the formula 



f dmexp(ikm(s — T))m a m b m z m r 
= ( s - ~ T )) _1 sin ( fc ( s " T ))- 


We have 
Hence, in eq.(35) 


Pjl;abPjl;ab — 6 . 


F(t,t';k) = IPji-jtirdtdt'fgdsf* ds'(a(t)a(t')) X Q 1 (k;t,s)G 1 (/c; ?, s')a{s)a(s') 
fg dr f Q s dr'Wji-j'i' (s — t,t,s' — t', t'; k)P G (k\ r, r'). 

(37) 

P(f , f ; k) is the spectral function for the magnetic contribution to the temper¬ 
ature fluctuations. For a small k it tends to g(t) (with a certain function g) 
because P a (k) —> Kf(s, s'), and the functions Wji-fi> and G~ 1 (k; t, s) in eq.(37) 
for a small k also tend to a constant multiplied by a function of time. For 
a large k the spectrum distribution P a behaves as P 7 (the one for gravitons). 

In the higher orders of the conventional perturbative calculations of the tem¬ 
perature fluctuations for a large k we would obtain (from eq.(19), with the 
energy-momentum T for matter fields on the rhs) the contribution to temper¬ 
ature fluctuations similar to the one resulting from eq.(37). Hence, we can 
conclude that the magnetic field does not substantially modify the behaviour of 
the spectral function for large k in comparison to the one without the magnetic 
held. Hence, it would not be detectable by a measurement of large multipoles. 

We do not have an explicit formula for Q^ 1 . However, for a small k , such that 
k « a~ 1 dfa, the dependence on k can be neglected. In such a case, eq. (37) 
gives an analytic formula for metric fluctuations caused by a linear dependence 
on the primordial magnetic held. After the analytic calculation of the spectral 
function F(t, t; k) in eq.(37) we are able to derive a formula for the contribution 
of the primordial magnetic held to the temperature fluctuations (26). We will 
be brief in the discussion of this derivation because it is already standard and 
described in many textbooks [20] [21]. We follow the calculations of our earlier 
paper [18] (concerning dissipative systems). In the integral dk = dkk 2 de we 
integrate hrst over e in the exponential in eq.(26). We obtain 


f de exp(— i(t — s)nk + i(t — s')n'k) 

= 2Trk~ 1 \(t — s)n — (t — s') n, | 1 sin (\(t ~ s)kn — (t — s')fcn'|^ . 
Next, we use the expansion 

fc _1 |(< — s)n — (t — s , ) n, | -1 sin (\(t — s)kn — (t — s')< 7 n'|^ 

= J2Zo( 21 + - s))ji(k(t - s'))Pi( nn'). 


(38) 


(39) 
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ji is the Bessel spherical function related to the Bessel function J [23] 


M z ) = (40) 

and Pi are the Legendre polynomials. 

If F (37) is known then owing to eqs.(26) and (34) there remains to perform 
the integrals over s and k 

(0(f, n)0(t, n')) = 

= i( 27 0“ 3 E“o fo ds fo ds ' f dkF ( s ’ s '> ^)(2(nA(k)n') 2 - (nA(k)n)(n' A(k)n')) 
ji(k(t - s))ji(k(t - s'))(2l + l)Pi(nn'), 

(41) 

where 

iiA(k)n / = nn' — k _2 (kn)(kn / ) = A(nn , 1 en, en), (42) 

nA(k)n = 1 — k~ 2 (kn) 2 = (5(en). (43) 

This formula is the starting point of calculations in [20] (see also our calculations 
in [18]). The expansion in Legendre polynomials reads 

(0(f, n)0(t, n')) = EEo( 2Z + l)Pi(t,nn')Pt(nn') = Y^Zoi 21 + l)C , ^(^)^ , /(nn , ), 

(44) 

(where Di is the term in front of (2Z + 1)P; in eq.(41)). In eq.(44) DiPi still must 
be expanded in Legendre polynomials if the coefficients Ci are to be independent 
of the angle. We have from eqs.(26),(36) and (41) 

A = fo d s fo ds' f 0 °° dkk 2 F(s, s', k ) 

(2A(nn', -id s ,id s ') 2 - d{-id s )S{id a ')^ji(k(t - s))ji(k(t - s'))- 

Let us consider only the term without derivatives in eq.(45) (denoted Di) re¬ 
sulting from the expansion 

2A(nn', —id Sl id s >) 2 - 5(-id s )5(id s >) = 2(nn') 2 - 1 + 0(d s ,d s >), 


where O is a polynomial of at least first order in derivatives. The terms in 
eq.(45) with derivatives can be calculated when Di are known [20]. We have 


Di = T A- T (2(nn') 2 - 1) f* ds' f* ds 
f 0 °° dqq 2 F(s, s', q)ji{q(t - s))j t (q(t - s')) 


(46) 


We make an approximation for F in eq.(37) which can be justified on the basis 
of the discussion following eq.(37) 


F(s,s',q) = g(s,s')q 3+e , 


(47) 


with a certain function g, where e is different for a large q and for a small q. 
The main contribution to the integrals (45) with the spherical Bessel functions 
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jl(rq) comes from rq ~ l. Hence, we can see that the behaviour of the integrals 
(45) corresponding to the spectral function F (47) at large q is responsible for 
the behaviour of Di at large l and the behaviour of F at small q corresponds 
to small l in Di. We can derive an exact result for the integral (46) using the 
formula 6.574 of [23] (7 < a) 


/ 0 °° d QQ 2 +€ Ji+i (eg) J; + 1 (79) 

F(l + f, + §,Z + §> ^2), 


1 _ ) 

2 ^r(|-f)r(;+|) 




(48) 


where F(a, 7 , z) denotes the hypergeometric function. Applying eq.(48) we 
obtain 


A(i) = ^(2(nn') 2 -1) 
g(s,s')F(l + f,-| + f ,1 


£d±j) 


r( |- §) rp+|) Jo ds /o ^ 

(t 


-s ; \ 


- s 


(49) 


The integral (46) can easily be calculated if F (47) is concentrated at s = 
s' = Sd- This case describes an instantaneous metric perturbation (the metric 
perturbation is limited to the moment Sd) corresponding to a sudden decoupling 
at s = Sd from the last scattering surface [20] [21]. In such a case s = s' = Sd in 
the argument of the hypergeometric function (49). We can obtain the value of 
the liypergeometric function at 1 using the formula 


F(cn, /?, 7 ,1) 


r(7)T(7-q-/l) 

r(7-a)r(7-/3)‘ 


Then, if e = 3 — p 

D l ~T{l+ e -){T{l + 2- e -))- 1 ~l 1 -o (50) 


Eq.(50) shows that the magnetic field is changing the behaviour of temperature 
fluctuations at small l. Di calculated with the magnetic field corrections do not 
decrease with l as they do for the inflationary P 7 ~ k~ 3 (then e = 0 or p = 3). It 
is known that there is a discrepancy between the theoretical and observational 
multipole contributions at low l to temperature fluctuations. A part of it can 
come from the primordial magnetic field. It follows that the part depending on 
the magnetic field behaves as 


(0(t,n)9(f,n')) B ~ ESo( 2l + “ ! ) 

^87rG(27r) _3 967ra(< ( i) _2 /3 _4 j b 2 (N + — AT) 2 d;Pilin'), 

with a certain slowly varying di. In eq.(51) td is the decoupling time, di can be 
calculated from eq.(37) and (41) by a numerical evaluation of the integrals. It 
varies slowly with l. There is still the contribution from the energy-momentum 
Tem °f the magnetic field (discussed in refs. [5]-[ 8 ]). This contribution depends 
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on b 4 G 2 a ~ 4 (the forth power of the magnetic field strength). It is decreasing 
faster with l because the power spectrum P 2 s(fc) is decreasing exponentially. 
The ratio of the numerical contributions of the linear and the quadratic ( Tem ) 
terms depends on several parameters : the total charge N + — N _, the strength 
of the magnetic field b and the temperature (/3a(t<j)) _1 at the decoupling. There 
are no precise estimates of these parameters. However, from the dependence of 
di on l we could infer the presence of the primordial magnetic field and the 
charge of the universe. 


5 Summary 

We have derived an elementary formula for the particle’s density distribution 
resulting from the perturbative solution of the Liouville-Vlasov equation. We 
have discussed a variation of the distribution which is linear in the magnetic 
field. It seems that this term has been ignored in the hitherto studies of the 
magnetic field in the universe. The linear term is non-zero if the primordial 
plasma is charged. There are strict estimates on the charge of the universe [9] 
[10]. The most elementary bound results from the argument that the electric 
repulsion cannot be much bigger than the gravitational attraction. This argu¬ 
ments restrict the ratio of the charge of the universe (in electron units) to the 
baryonic number to be of the order 10~ 18 . The strength of the magnetic held 
is also restricted to be extremely small: 1.0 nG in the epoch of the photon last 
scattering [1], As a consequence the linear term gives a small contribution to 
the temperature fluctuation spectrum. We have calculated its dependence on 
the magnetic held. There are some undetermined parameters in the formula. 
However, the functional form of the temperature fluctuations could discriminate 
between various models of the primordial magnetic held. The usually discussed 
quadratic perturbation of Einstein equations resulting from the electromagnetic 
energy-momentum has a contribution to the temperature fluctuations which de¬ 
pends in a different way on the probability distribution of the magnetic held and 
it does not depend on the temperature at the decoupling. 
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6 Appendix 

We wish to calculate 


A(k) = J dqP 7 (k - q)P B (q) 


( 1 ) 


with 

P 1 {k) = k~ a 


( 2 ) 
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as a typical model for a graviton distribution and 


P B (k) = k a ex p(-y) (3) 

for the magnetic field probability distribution. Then, eq.(l) gives 

Pcr(k) = 2ir / Q °° dq dxq s+2 (k 2 + q 2 — 2qkx)~% exp(—^-) 

= 2irk~ 1 (2 — a) -1 / 0 °° dqq s+1 + q) 2 ~ a — \k — q\ 2 ~°^j exp(—^-) 

= 2nk~ 1 {2 — a) -1 / Q fc dqq s+1 ^(k + q) 2 ~ a — (k — g) 2_a j exp(—^-) (4) 

+2irk~ 1 (2 — a) -1 f£° dqq s+1 k + q) 2 ~ a — (q — k) 2 ~ a J exp(—^-) 

= J 1 (fc) + J 2 (fc). 

It is easy to see that 

h(k) = k~ a g(k) (5) 

where limfc_>. 00 g(k) = const ^ 0 and I 2 decays exponentially for a large k. 
Moreover, for a small k we obtain that P CT (0) > 0 is finite if 5 — a > —3 . 

Concerning the convolution (34) of the magnetic spectral function it is useful 
to consider its Fourier transform P B { y) 

Pb{ k) = J dyP B {y) exp(iky). (6) 

Applying the Fourier transform for the distribution (3) with o = 0 we obtain 

P B {k) = exp(-y) = J dyexp(-^y 2 )exp(zky)(7rA)"2. (7) 

Hence, the spectral function defined in eq.(34) is 

P 2 B(k) = J dpP B {p)P B (k — p) = (|)" 3 /dyexp(-|y 2 )exp(tky) 

J dy'6(y — y') exp(—^y' 2 ) ex p(iky / ) = 8(27 tA) 2 exp(— 1 j 

For general er yk 0 it is difficult to obtain explicit formulas for P 2B . In general, 
in eq.(6) if P B is decaying faster than any polynomial then also P B is decaying 
in the same way and the convolution (34) of such functions has this property. 
From eq.(6) we can also derive the relation 

P 2 B(k = 0) = |dy(P B (y)) 2 . (9) 

Hence, if P B ( y) is square integrable then P 2 ^(k) is regular at k = 0 and disap¬ 
pears fast at large k. The behaviour for small k is not generic with respect of 
convolutions and Fourier transforms. We do not have special reasons to apply 
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the spectral functions exactly of the form (3). In the literature a sharp cut-off 
at y/X is also applied [24]. For the properties of the stochastic magnetic fields 
only the behaviour of Pb for a large k and a small k is essential. We suggest to 
replace the spectral function (3) by 

/ o° 7.2 

exp(—s—). (10) 

This function has the same behaviour as the one in eq.(3) for large as well as 
small k. It has the virtue that its Fourier transform can easily be calculated 
and we obtain a workable representation of P‘2B ■ It can be shown by an explicit 
calculation that if a > —1 then P2s(0) ^ 0 is finite . It decays exponentially for 
a large k. If a < —1 then P 2 s(fc) is singular when k —> 0 but its singularity is 
less than k~^ ( e.g., as discussed in [24], when a = — 2 then we have P 2 B — k~ l 
for a small k). 
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